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We present a phenomenological study of Drell-Yan pair production at hadron colliders based 
on the NNLO fixed order calculation and on NNLL resummation of threshold logarithms. We 
give an argument to prove that resummation effects are relevant also for values of x = M^/s far 
from threshold. We compare different prescriptions for the calculation of resummed quantities, 
emphasizing the differences coming from subleading terms, which are important when x is small. 
We present phenomenological predictions for Drell-Yan rapidity distributions at the LHC, we 
study the ambiguity related to the resummation prescription, and we compare it to that coming 
from scale variation. 
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1. Introduction: threshold resummation in the Drell-Yan process 

A generic (differential) parton model cross section a at hadron colliders can be written as a 
convolution 

dz ^ f 



where x = M'^ /s (M is the Drell-Yan pair invariant mass, ^/s the hadronic cms energy), L(z) is a 
parton luminosity and d(z) is the parton-level cross-section. In d(z) logarithms of 1 — z appear at 
all orders in ttg 

"log'(l-z)" 



a: , I = 2n-l,...,0 , (1.2) 

+ 

and in the partonic threshold limit z — )• 1 these logarithms become large, spoiling the reliability of 
the perturbative expansion. All order resummation [1] is needed in this case. 

2. Relevance of resummation of log( 1 — z) at small x 

We see from eq. (1.1) that the partonic threshold region z — )• 1 is always included in the convo- 
lution integral. To understand when resummation is relevant at the hadron level we need to establish 
when z — )• 1 region gives the dominant contribution to the integral. Qualitatively [2] this happens 
when the shape of the parton luminosity enhances that region of integration. A quantitative answer 
[3] can be given working in N -space, where N is the Mellin-conjugate variable to x. Eq. (1.1) can 
be rewritten as a Mellin inversion integral 

a{x) = ^i dNx-^L{N)a{N) = :^. emogi+iogL(N)+ioga{N) ^2.1) 

where L{N), d{N) we. the Mellin transforms of L(z), o{z). This inversion integral is dominated 
by the region where the exponent has a minimum in the positive real axis (saddle point A'^q). By 
a general argument, one can show [3] that a saddle point always exists. We show in Fig. 1 the 
position of the saddle point A'^o as a function of x, where d{N) is the DY qq contribution at NLO 
and the paiton luminosity is computed using CTEQ6.6 in pp collisions. The saddle point A^o is 
monotonically increasing with x, meaning that at large x the contribution to the cross section mainly 
comes from the large A'^ region, which corresponds to the partonic threshold region. For smaller 




Figure 1: Position of the saddle point A^o a function of x for the order-a^ Drell-Yan qq cross section. 
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X, the saddle point Nq decreases; to evaluate when the large logarithms still give the dominant 
contribution, we use a fixed order computation to compare the logarithmic terms with the full 
result. The DY qq contribution at the NLO is given by (up to an electroweak normalization factor) 

a(z) = 5(l-z) + ^Ci(z) + 0(a2) (2.2) 

In Fig. 2 we show the Mellin transform of C\ (z) (black solid curve) and its logarithmic part, the 
first term in eq. (2.3) (blue dotted curve). The red dashed curve corresponds to the first two terms 
in eq. (2.3): indeed, the term has the same dynamical origin of the logarithm inside the plus 
distribution, and may be included as well in resummation (see Sect. 3). Inspection of Fig. 2 leads to 
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Figure 2: NLO Drell-Yan coefficient function as a function of A^, and its logarithmic approximations. 



the conclusion that, down to values of N around 2, the logarithmic part dominates the partonic NLO 
cross section. Hence, if the saddle point A'o > 2, the integrand defining the hadronic cross section 
o{x) in eq. (2.1) is dominated by the logarithmic part of the partonic cross section. Comparing 
with Fig. 1 we can conclude that for values of x > 0.003 resummation of the logarithms is relevant, 
in the sense that it includes to all orders the dominant terms of the perturbative coefficients. Note 
that this value is rather small, and usually considered too small for threshold resummation to have 
a sizable effect. 

3. Comparison between the Borel prescription and the minimal prescription 

Threshold resummation is performed in A'^ space [1], and the generic resummed quantity can 
be written as a function of log ^ {h^ are coefficients which depend on «§) 

^res = £ log - j = £ oc' log' - , « = 2j3oas • (3.1) 

This series has a radius of convergence 1 in a log ^, since the sum has a branch cut in <N < +oo 
due to the Landau pole of OCg- For this reason, the Mellin inversion integral (2.1) is not defined, 
because the parameter c should be to the right of all the singularities of the integrand, but the cut 
makes this condition impossible to be satisfied. Otherwise stated, the series obtained performing 
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the Mellin inversion of (3.1) term by term is divergent [4, 3]. Hence a prescription is needed to 
obtain a resummed cross section in the physical x space. 

In [5] the minimal prescription has been proposed: the idea is simply to choose the integral 
path in (2. 1) to the left of the cut. This choice has some good properties, one of them being that the 
result is asymptotic to the divergent series. However, it also has some undesired features: first it is 
defined only at the hadron level, and the integral gets a contribution from an unphysical region of 
the parton densities [5]; second, there is no way to control subleading terms. To clarify this second 
statement, we apply the minimal prescription to the ^-th term of the series (3.1), i.e. log*^ ^. In this 
case the integral is an exact Mellin inversion, and the result is a sum of terms of the form 



log^logi 



(3.2) 

1. There is no way to 



which are an approximation of the physical logarithms (1.2) in the limit z 
force the minimal prescription to reproduce the physical logarithms. 

In [4] a different prescription, based on the Borel summation of the divergent series (3.1), has 
been proposed. We present here a simpler proof of the Borel prescription formula, following [3]. 



First, we rewrite the Mellin inversion of log*^ as 



N 



k\ 
2ni 



log 



(3.3) 



where the integral path in the r.h.s. is any closed contour which contains i§ = 0. Next, we use the 
identity 

k\= r dwe-'^'w^ (3.4) 



to eliminate the k\ in (3.3). The Borel method consists in exchanging the w integral with the series 
in (3.1), obtaining 



1 

2ni 



d^ 



+=" dw 
a 



-w/a £ 



w 



log^-^ 



1 



(3.5) 



However, this expression is still divergent. The Borel prescription is obtained putting an upper 
cutoff C to the w integral; this choice amounts to the inclusion of an higher twist term. It is shown 
in [4] that this result is asymptotic to the divergent series. Most importantly, the z dependence is 
explicit, and can be modified in order to better reproduce the physical logarithms. Since in the 
previous section we have seen that the inclusion of the log ^/z. term can help to better approximate 
the shape of the full fixed order, we include also these logs and obtain our final expression: 



1 



d^ fCdw 



iTiiJ t,T{^) Jo a 



M'/a ^ 



w 



(3.6) 



The main features of the Borel prescription are that it is a paiton level expression and that we can 
control the z dependence. Furthermore, the parameter C can be used to estimate the ambiguity of 
the prescription. 
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Figure 3: Rapidity distribution of neutral Drell-Yan pairs of invariant mass M — mz (upper plots) and 
M = 1 TeV (lower plots) in pp collisions at ^ = 7 TeV (using CTEQ6.6 pdf set and a^{mz) — 0. 118). 



4. Results: Drell-Yan rapidity distributions at hadron colliders 

In Fig. 3 we show the theoretical predictions for the production of a neutral DY pair in pp 
collisions at y^s = 7 TeV. In all plots, the 7 < region shows the fixed order results, at LO, NLO and 
NNLO accuracy. The Y >Q region includes the resummed results (using C = 2 Borel prescription 
in the left plots and minimal prescription in the right plots) at LL, NLL, NNLL accuracy. The error 
bands are obtained by varying the renormalization and factorization scales independently by a fator 
of 2 around M. The pdf uncertainty is not included. 

The upper plots refer to the production of a DY pair of invariant mass M = mz- In this case 
X ~ 10^"^ is rather small (with respect to the value 0.003 found in Sect. 2), meaning that the large z 
region of the partonic cross section does not give the dominant contribution to the cross section. In- 
deed, there is a quite large difference between the Borel prescription and the minimal prescription 
results, coming from the different way of including subleading terms. We believe that this dis- 
crepancy can be used in order to better estimate the ambiguity due to the unknown higher orders. 
Moreover, the resummed contribution does not reduce the scale dependence, since resummation 
only affects the qq channel, and the reduction of scale dependence is a combined effect of the qq, 
qg and gg contributions. 

The lower plots refer to the production of a DY pair of invariant mass M = I TeV. In this 
case, X ~ 0.02 is in the region for which resummation is relevant. At the NNLO-i-NNLL level of 
accuracy, the inclusion of the resummed term gives a small improvement: the scale uncertainty is 
mildly reduced. 

5. Conclusions 

In conclusion, we have shown by a quantitative argument that inclusion of threshold logarithms 
at all orders does improve the accuracy of QCD perturbative predictions even when x = M^/i is not 
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close to 1 (down to x ~ 10 ^ — 10 ^). Furthermore, we have presented a realistic phenomenological 
application of the Borel prescription for the computation of resummed cross sections. 
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